THE RADIATIVE DECAY OF VECTOR MESONS 

Ting-Liang Zhuang 
Center for Nonlinear Science, University of Science and Technology of China 
Hefei, Anhui 230026, People's Republic of China 

Xiao-Jun Wang 

Center for Fundamental Physics, University of Science and Technology of China 
Hefei, Anhui 230026, People's Republic of China 

Mu-Lin Yan 
Center for Advanced Study 
Tsinghua University, Beijing, 100084, People's Republic of China 

and 

Center for Fundamental Physics, University of Science and Technology of China 
Hefei, Anhui 230026, People's Republic of China]] 



In this paper, radiative decays p° — > 7r + 7r~7, 7r°7r°7 ,0 — > K + K~ 7 , K° 'if 7 are studied 
systematically in the U(3)_l xU(3)_r chiral theory of mesons. The theoretical differential 
spectrum with respect to photon energy and branch ratio for p° — > 7r + 7r _ 7 agree well with 
the experimental data. Differential spectrums and branch ratios for p° — > 7r°7r°7, (f> — > 
K+K--y,<f) -v if°i?°7 are predicted. The process <f> — > K°K°j is relevant to precision 
measurment of CP-violation parameters in the kaon systerm at a (^-factory. We give a 
complete estimate of the branch ratio for this decay process by including scalar resonance 
/o,ao poles, nonresonant smooth amplitude and an abnormal parity process with K* pole 
which hasn't been considered before. We conclude that processes with intermediate K* do 
not pose a potential background problem for <f> — > K°K° CP violation experiments. 
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I. INTRODUCTION 



Radiative decays V — ► PP'y(wheie V denotes vector mesons, P denotes pseudoscalar mesons) have at- 
tracted much interest in past decade[2-10]. The study of this kind of rare decays is important in hadron 
physics, both because it is intimately related to the QCD-inspired descriptions of the_dynamics of mesons 
and because it has been urged by experiments. For example, the reaction <f> — > K°K°j poses a possible 
backgroud problem of <p — > K°K Q at future cf> factory. The latter process has been proposed as a way to 
study CP violation JjJ. 

The purpose of our present paper is to systematically study the processes p° — * 7r + 7r~7, 7r°7r°7, 4> — > 
K + K~~/, K°K°7 in the framework of U(3)z,xU(3).r chiral theory of mesons E2|. In this effective chiral 
model, all couplings among pseudoscalars, and its lowest resonances are fixed by introducing an universal 
coupling constant g. Thus, a unified description of mesons in the low energy is provided. 

In fact, this effective model is an extended chiral quark model including — , 1^ mesons. The chiral quark 
model, originated by Weinberg and then developed by Manohar and Georgi fTi) ], provides a QCD- 
inspired description on the simple constituent quark model. In the view of Manohar- Georgi model, between 
the chiral symmetry breaking scale(A xSB ~ 1 — 2GeV) and the confinement scale {h-QCD ~ 0.1 — 0.3GeV), 
the dynamical field freedom are constituent quarks (quasi-particle of quarks), gluons and Goldstone bosons 
associated with chiral symmetry spontaneously breaking. In this quasiparticle description, the effective gluon 
coupling is small and interactions between quarks and Goldstone bosons is important. The external gauge 
fields(e.g., photon field) can be introduced by localizing the global chiral symmetry. On the other hand, it 
is well known that in the electromagnetic interaction of mesons, the vector mesons play an essential role 
through VMD (Vector Meson Dominate) fllSj , Therefore, it is quite nature to extend chiral quark model to 
include spin-1 meson resonances via VMD and via minimal coupling principle. 

The U(3)z,xU(3)fl. chiral theory of mesons has been studied extensively [|ll||l5| , |l6f| . The basic inputs for 
it are the pseudoscalar decay constants fp, vector meson mass my and a universal coupling constant g. 
Predictions of this model are in good agreement with data |lj,|l5],|l6| . In particular, in Ref. [jljj it has 
been shown that the low energy limit of this theory is equivalent to the chiral perturbation theory, and 
the QCD constraints in Ref. |19) are satisfied by this model. Therefore, as an effective model of QCD, the 
U(3)l xU(3)i? chiral theory of mesons is reliable. 

The content of this paper is organized as follows. In Sec. 2, we present a brief review of chiral quark model 
and the basic notations of the U(3)lxU(3)h chiral theory of mesons. In Sec. 3 and Sec. 4, the branch ratio 
for these decays are calculated respectively, and the gauge invariance of these decay amplitudes is checked 
explicitly We give a summary of the results in Sec. 5 



II. CHIRAL QUARK MODEL AND U(3) L xU(3) H CHIRAL THEORY OF MESONS 

The simplest parametrization of chiral quark model is 

C x = tp(x)(ij ■ d — mu(x))ip(x) (1) 

where 

tp = (u,d,s) T , 
u{x)= l -{l- lb )U{x)+ l -(l + 1 s)U\ 
U{x) = exp(2a a $ a // ), 

and A a are Gell-Mann matrices of SU(3), <& a are fields of pseudoscalar meson octet and /o ~ / w = 186MeV, 
m is a parameter related to the quark condensate. 

This Lagrangian is invariant under global chiral symmetry transformation. The vector(V M ) and axial- 
vector^^) external field are introduced into C x due to the requirment of local chiral symmetry, i.e., we can 
replace the derivative operator in Eq.(l) by covariant derivative operator with affine connectionfor gauge 
potential) V p + .4^75 as follows: 

d„ -» = d„ ~ i(V„ + Atfs). (2) 
Then we obtain a theory describing strong and elecro-weak interactions of mesons. 
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Spin-1 meson resonances can be included via VMD, i.e., via substitution of a new affine connection 
((V + u) M + (A + a) M 7 5 ) for the former one V M + A^/ 5 in Eq.(2), 

W V M = - i((V + «)„ + + a) M75 ), (3) 

with 

2 1 11 

-^As)/, + (- - 

,2 1 . ,1 1 

3 + Vl A8K + ( 3 " W 



V = Ti< + A a iff + (- + — X S )U + ( o ~ — -W^, 



7*pj, + A Q if; Q + (- + — A 8 K + (- - — A 8 )<^, (4) 



where i=l, 2, 3 and a=4, 5, 6, 7. 

Now, the author of Ref. jl2j came to this extention and proposed a Lagrangian, 

C = i)(x)(i~f ■ (d - i((e QA + v) + aj 5 )) - mu(x))ip(x) + ^ml(v^ a v% + a M0 a£), (5) 

where is photon field, Q is the electric charge. The mass term of v and a is chiral gauge invariant because 
the v and a transform homogeneously under local U(3)l x U(3)r symmetry. Note that there are no kinetic 
terms in Eq.(5) for all meson fields, since they are treated as composited fields of quark fields instead of the 
fundamental fields. The kinetic terms for these fields will be generated via loop effects of quarks. 

Following Ref. p2|, the effective Lagrangian of mesons (indicated by M) are obtained through integrating 
over the quark fields, 

exp{i J d 4 xC M } = J [#][#]exp{i J d 4 xC}. (6) 

Using the dimensional regularization, and in the chiral limit, the effective Lagrangian Cre (normal parity 
part) and Cjm (abnormal parity part) have been evaluated in Rcfs. [ fL2| . The Lagrangian describing normal 
parity processes reads 

Cre = ^TrD^UD^ - ^Tr(L^L^ + R^R^) 



+i^^Tr{D lx UD v U^L Vil + DpDvUR^) + ^^TrD^UD^D^ 

N c 



T Tr{D^UD^D v UD u U^ + D^UD U U ] D V U - D fl UD v \J^ D^UD^U^) 



12(4tt) 2 " 

■~m 2 Tr(L^ + R^), (7) 



where 





= - a 










- iL^U + iUR^, 






= v^ f 








= S7/j,L v 


~~ V u-Lfj, — i[Lfj,, L v \ -f- 






= S7nRv 


— SJpR^ — i[Rp, Rv\ - 


f eoQF* , 




= d^ + 


[—ieoQAft, # = 


u,u\l,r, 


F L 


- F R - 


d^A v — d v A^ . 





(8) 

Here an universal coupling constant g has been introduced to absorb the logarithmic divergence due to the 
integral of quark loop. 

In Lagrangian ( |?]) the field a^(x) mixes with d^lx), which should be diagonalized conveniently via field 
redefinition, 

g 

R^ R^ + i-U^^U . (9) 
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Then the following equations are derived, 



F 2 2c f 2 

) = 1, c = for two-flavor. (10) 

f4 9 2gm P 

F 2 2c' f 2 

■75- (1 ) = f , c' = ^ , for three-flavor (f 1) 

It should be stressed that this field redefinition is different from the one in Ref. |fl^] , Op — > a M — f^^- 
Eq.( ^) keeps the chiral symmetry explicitly. It plays an important role when we check the electromagnetic 
gauge invariance of the decay amplitude in following sections. Eq.(7)-Eq.(ll) provide the formalism employed 
in this paper and all the calculations are performed in the chiral limit. 



III. THE DECAY p° -> tt+tt^, 7r°7r°7 

In this section, we will restrict our calculations to the two-flavor case. The needed vertices to evaluate 
these processes can be obtained from Sec. 2. For the reaction p° — > 7r + 7r~7, the Feynman diagrams are shown 
in Fig.nl the involved vertices are: 




FIG. 1. p° 



IT IT 7 



where 



C p7I ^ = eA p p„AW + eA^d^F^wV , (12) 

Cpx* = AfryWey* , (13) 

+B p P y v d^ + B^p k XV], (14) 

£ 77rT = eA^d^^e^k , (15) 

= e ^-^=^ Fu, u a^d"Tr j e i jk , (16) 
2gTT z b z 



9 l 2** ft" 9' 2 
B? = - A s + -g^, oP_j^^ 3/, 



RP = J 2 ^ RP = i 

7=(l-m) _i ( 17 ) 
The decay amplitude of this process is: 

M 

= M a + M b + M c , (18) 

where 
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M a = < n + n-j\iT J d 4 x£ p7I ^{x)\p° >, (19) 

M b = < 7r+7r-7|t 2 T J d*x J dtyCp^toC^Wlp >, (20) 

M c = < 7r+7r- 7 |z 2 T J d 4 x J d 4 yC pa7T {x)C^{y)\p° > . (21) 

Since photon is on-shell in this process, we should check the gauge invariance of this amplitude. From Eq.(16) 
we see that the vertex of C ia -R is already gauge invariant, therefore, we need to check only the sum of M a 
and Mb- One can derive from Eqs. (19-20) that 

M a = 2ie(A{g liV + A p 2 p ■ qg^ - A^p^s^e- , 

Mb = 2ieAl{^ + , (22) 



where p,k + ,k ,q denote the momenta of p°, ir + , it and photon fields respectively, and e^, e~ are the po- 
larization vectors for p meson and photon field respectively. If we substitute with q v in Eq.(22), we will 
obtain: 

{M a + M b )\e^ r cc (q, + k+ + k-)e$ . (23) 
Using four-momentum conservation and the space-like condition of the wave function of vector field, 

p = q + k + + k~ , 

V ■ ef = • (24) 

Eq.(22) vanishes, so the gauge invariance of this amplitude is kept. 

Before we give the numerical results of the width and the branch ratio of this process, it is necessary 
to point out that there is no adjustable parameter in our calculation. The basic input g, as an universal 
coupling constant in this theory, can be fixed by a experiment, thus we can compare our theoretical results 
with experimental data and give predictions for processes which have not been measured in experiments. In 
this paper, we take g = 0.39. Then we obtain, 

r(p° -► TT+TT^) = 1.54 MeV , 

B(p° -> 7T+7r"7) = 1.03 x 10~ 2 for E 1 > 50 MeV. 

which compares favourably with the experimental data f^O), B exp (p a — > 7r + 7r~7) = (0.99 ± 0.04 ± 0.15)10~ 2 
for E 1 > 50MeV, where E 1 is the photon energy in the rest frame of p° . The shape of differential spectrum 
with respect to photon energe is also given in the Fig. ^[ One can see that the experimental result is in good 
agreement with our theoretical expectations. 

The reaction p° — > 7r°7r°7 involves only abnormal parity, the Feynman diagrams is shown in Fig.^|. Fol- 
lowing Ref. fll2||, by means of bosonization the quark propagator, we obtain 




FIG. 2. p° -» 7r°7r°7 

/W = Y^e^d^p^dpir 1 , (25) 

9 Jit 

C-t^ = ~e-^—e^d^ v A a dp^ . (26) 
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The vertex £ 7 ^ w is obviously gauge invariant because of the totally antisymmetry tensor £^"-0 . The branch 
ratio is calculated to be 



B( P " -> 7T U 7T U 7) = 1.02 x 10" 



(27) 



and the differential spectrum with respect to photon energy is shown in Fig.|^ which will be tested in future 
experiments. 



IV. THE DECAY 



KK 1 



In order to calculate the <fi decay, we need to extend our calculation to three-flavor case. First, we derive 
the vertices for the transition <fi — + K + K~"f, the Feynman digrams are shown in Fig.[|, the vertices of this 
process have both normal parity part and abnormal parity part. The normal parity part is 



where 



C 



<pK+K~j 



A^K+K- 



C*k+k- = iAf^K+d^K- , 
£ yK +K- ^leA^K+d^K- , 

iBi<t>^ v d^K + + iBid 2 ^KrK+] + h. 



c 



~/KiK 



3 

2gir z r z 



V2ml 



f 2 



and the abnormal parity part is 

£<pK*K± 



r' 1 1c' 

[4c'(l--) + ^(l--) 2 ], 
k 9 9^ 9 

y^B?(c^c',U^ f K ), I = 1,2,3, 4. 



3y/2 



^d u 4> a K*+d K- + h.c. , 



e^ a(3 d v A a K* + d p K- + h.c. 



(28) 
(29) 
(30) 

(31) 
(32) 



(33) 

(34) 
(35) 





FIG. 3. 



K + K~ 



The gauge invariance can also be checked as in the p decay. Using the same value of g, i.e., g = 0.39, we 
give the differential spectrum in Fig.|| . 

Now, we consider the transition <fi — > K K°-f. The photon in this reaction is very soft(E max < 25Mev), 
so it is difficult to distinguish it from the genuine <j> — > K°K° events which has been proposed as a way 
to measure the small parameter in studying CP violation. The branch ratio of <j> — > K°K°j(^ 10~ 6 ) will 
limit the presion of this measurment. This quantity has been predicted by several authors [4- 10]. In Refs.[4- 
8], the contribution of interchanging the scalar meson S(o,q,ot Jo) has been obtained via chain reaction 
4> — > S*7 — > K°K°j, in which the decay <j> — ► S7 is proceeds through the charged K loop. The uncertainty of 
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this approach is that the coupling constant ggxR 1S n °t wen know because of the lack of experiment data. 
In Ref. [9fl, the non-resonant contribution has been calculated using current algebra and low energy theorem. 
In Ref. [1(J, the authors did not introduce ao, /o explicitly, they calculated the final state intereaction of KK 
system in a chiral unitary approach. This approach generates the ao,/o meson dynamically, the obtained 
amplitude after summing over an infinite series diagram also contain non-resonant contribution. All these 
calculations, however, did not contain the contribution of an abnormal parity process via interchanging K* 
which in principle must be added to the resonant poles and the nonresonant smooth amplitude, otherwise one 
may not assume scalar meson dominance a priori. In the present paper, we provide a complete calculation 
on the branch ratio for this decay by including the abnormal parity process with K* pole, scalar resonance 
/o, do poles and non-resonant amplitude. The role of scalar resonance will be dealt with as in the former 
works [4-8, 10]. The difference between their scheme and ours is that the needed vertices to calculate the 
loop diagram </> — > K + K~ r y has been derived in Eqs. (28-30) as well as all the coupling constants has been 
fixed by the universal constant g, in other words, there is no adjustable parameter in our calculation. The 
complete intereaction of this process including normal and abnormal parity vertices is: 

C = eAf^A^K+R- + iAf^K+d^R- + ieA^K+d^R- + L^k-k + £^k*k° (36) 

where 

2^ 2 J K 



£4>k*k° = - 2 2 f e^d u <t> a K;%K° + h.c. , (37) 



£ iK *ko = -^^e^d v A a Kfd (j K* + h.c. , (38) 

described the abnormal parity process. Denoting the momenta of K°,K° as ki,k 2 and defining si = 
(ki + q) 2 , s 2 = (k 2 + g) 2 , we derived the amplitude for this abnormal diagram(shown in Fig.||): 




with 



FIG. 4. (t> -> K°K°j through K* . 



3\/2 

^ 9 J K 



3\/2" 

ie J W72 ( d og^ + dikipkiv + dikz^kzv + d 2 k ltl k 2v + d 2 k 2ll k lv ) eSe£ (39) 



d ° = 77 2~r( 4 ^ 1 'hs + ai — mi) + (ki <-> fc 2 ), 



S _l (Al. I. I -' 

4(s 1 -mf f .) 
1 , s 2 ml- si 



di = -( 5 - — =— ), di = di(ki f-> k 2 ), 

I s 2 — rn K , s i — m K* 

d 2 = ^—r 1 o , d~2 = d 2 (ki <-> k 2 ). (40) 

s 2 - m K , s 2 - m K » 

In order to derive the contribution of scalar resonance /o and do, we need to calculate the one loop diagram 
though K + K~ and the final state intereactions of K + K~ to K°K° shown in Fig.|| .Similar to Ref. jl(|, we 
get the amplitude: 
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/ K + K .Jl^ K' I 

v -'-ex kx 



K" 

(a) 





\ 



(b) 

FIG. 5. (a).</> — > K°K°-y through charged K loop. (b).A" + _R'~ — > A' ,/? amplitude, the intermediate loops contain 
7T7T, K K, ir°rj. 

M ? = o 2 ' 1 2 7 ( fl ' & )b ' ?#/^ - P*^] £ ff e f ** ( 41 ) 

where a = m|/m|., 6 = Q 2 jm 2 Kl Q 2 = {k\ + k 2 ) 2 , 

I(a,b) = - J— - - /(I)) + -^(5^) -<?(")) (42) 

2(a — oj (a — o) 2 o a (a — o) z b a 



with 



{— arcsin(^^) 2 x> j 



_ J (4a:- 1)5 arcsin^) 2 a; > | 
9[X) -\ 4# [ln(2±) 



1 



J?± = — (l±vT^4i) (43) 



the t s (Eq.(9) in |10|) is the scattering amplitude of K + K~ to K°K a , its diagrammatic meaning is shown in 
Fig.^(b), and its explicit expression can be obtained from Eq.(30) in Ref. [ p"l| which is derived from Lippmann- 
Schwinger equation in the coupled channel approach. As declared in Ref. |Q, this amplitude(Eq.(41)) also 
take into account nonresonant contribution. 
The width and the branch ratio are given by: 

T(<p -> K°K^) = — - 2 3 / d Sl dQ 2 (M i + M 2 ) 2 (44) 
1927r^m^ 7 

B((j> -> K°K Q -/) = T(0 -> K°K -/) /4.43 (45) 

where a = e 2 /47r = 1/137. In the following numerical evaluation, the constant <? take the same value 
0.39 as in the previous. If we neglect the abnormal parity process, i.e. set Aii = 0, then we obtained 
£?(<ft — > lf i? 7)«caZar = 5.6 x 10~ 8 , which is a little different with the value 5 x 10 -8 quoted in Ref. 
pffl. On the other hand, if we neglect the contribution of scalar resonance, i.e. set M. 2 = 0, we obtained 
B{<f) — > K°K° 7) abnormal = 7.6 x 10 -8 . We see that the contribution of this abnormal parity process is the 
same important as the scalar resonance poles, so its contribution can not be neglected. After performing the 
integral of Eq.(44), we obtained: 

B{(/) -> K°K°j) = 1.8 x 1(T 7 . (46) 

We see although the interference is constructive, it will not provide much significant background for precision 
test of CP-violation in <f> — > KK. The differential spectrum with respect to photon energy for these three 
case(abnormal, scalar, interference) are given in FigM 



8 



V. SUMMARY 



To conclude, in this paper, we perform a systematic calculation of p and fts radiative decays in an extended 
chiral quark model, in which all coupling arc fixed by the universal coupling constant g. The gauge invariance 
of these decay amplitudes has been checked. The theoretical differential spectrum with respect to photon 
energy and branch ratio of p° — > tt + -k^j agree with the experimental data well. Predictions of differential 
spectrum and branch ratio_for the processes p° — > 7r°7r°7, <j) — * K + K~j, K°K°j have been derived. The 
branch ratio for <f> — > K°K°~f including the contribution of abnormal parity process with K* pole, scalar 
resonance ao,/o poles and nonresonant amplitude has been calculated to be about 1CP 7 and will not limit 
the precision measurment of the small CP-violation parameters at future cf) factory. 
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FIG. 8. Photon spectrum, dF/dE-,, for the process <f> — > K + K~>y. For > 5MeV 
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FIG. 9. Photon spectrum, dF/dE^, for the process -> K°K a j. Dot line: distribution only taking into account 
contribution of abnormal parity process with K* poles, dash line:distribution only taking into account the contribution 
of scalar resonance poles and nonresonant amplitude, solid line: the total distribution. 
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